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Abstract 

Cytoplasmic incompatibility (Cl)-inducing endosymbiotic bacteria, such as Wolbachia and Cardinium, have been well studied 
through field data and validations on the basis of numerical simulations. However, the analytically derived equilibrium 
frequency of multiple infections has not yet been determined, although the equilibrium for cases of single infection has 
been reported. In this study, we considered the difference equation for endosymbionts using three parameters: the 
probability of the failure of vertical transmission {/<), CI strength (-), and the level of host inbreeding (/>). To analyze this 
model, we particularly focused on Q^v/ i-^., the frequency of host individuals completely infected with all A^-bacterial strains 
in the population. Ql/, Qn at the equilibrium state, was analytically calculated in the cases where N=l,2 and is any 
arbitrary value. We found that can be described using two parameters: N and a, which is identical to i-i/pz. Ql/ has a 
larger value in a system with a smaller a. In addition, a determines the maximum number of strains that infect a single host. 
Our results revealed the following: i) three parameters can be reduced to a single parameter, i.e., a and ii) the threshold of 
the maximum number of infections is defined by a, which prevents additional invasions by endosymbionts. 
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Introduction 

Endosymbiotic bacteria such as Wolbachia and Cardinium are well 
known to be reproductive manipulators infceting insect cells [1,2]. 
Cytoplasmic incompatibility (CI) is considered to be the most 
common and efficient form of manipulation that can spread an 
infection throughout the host population. CI causes males infected 
with a bacterial strain reproductively incompatible with uninfected 
females [3]. Because of CI, infected females exhibit a relatively 
higher fitness than uninfected females; therefore, the number of 
infected individuals gradually increases and they become domi- 
nant in the host population. For example, approximately only 
10% of Drosophila simutans individuals were infected with Wolbachia 
in the middle of the 1980s in California; however, the infection 
rate increased to 95% in 1993 [4,5]. In a rearing system of Encarsia 
pergandiella infected with Cardinium, the infection gradually spread 
to almost fixation within the population, irrespective of the initial 
rate of infection [6] . 

The dynamics of Cl-inducing bacteria can be primarily 
determined by three parameters: vertical transmission efficiency, 
fecundity of infected females, and CI levels [7], although the 
dynamics may be affected by other parameters [8,9]. Vertical 
transmission, from a mother to her offspring, is the main route of 
transmission of the bacteria to other host individuals. The 
transmission rate is defined as the proportion of offspring infected 



from an infected mother, which is approximately 1 (> 0.95 in most 
cases) [10]. Successful transmission has crucial effects on the 
spread and maintenance of the bacteria [1 1]. Cl-inducing bacteria 
sometimes affect the fitness of the infected host diretctly through 
fecundity. The direction of their effect through fecundity depends 
on the host-bacteria combination, and the effect may be negative 
[12] or positive [13], or often neutral [14]. The CI level is defined 
as the proportion of the number of offspring died because of CI 
relative to the number of total offspring reproduced. In addition, 
the effect of CI is dependent on the host-bacteria combination, 
which highly varies from zero to one, i.e., from nearly neutral to 
complete mortality [15-17]. 

Models have been developed to study the dynamics of Cl- 
inducing bacterial symbionts, particularly in Wolbachia, using the 
three conventional parameters [4,18]. The long-term behavior 
depends on the initial frequency of bacterial infection [19]. If the 
initial frequency of bacterial infection is below the threshold, the 
frequency heads towards extinction; however, if the frequency 
stochastically exceeds the threshold (e.g., because of random 
genetic drift), then it is expected to spread to another equilibrium 
state where both infected and uninfected individuals exist. A 
higher vertical transmission rate and/or stronger CI will lead to 
higher infection frequency equilibrium. In most of the models on 
the basis of difference equations, the equilibria have rarely been 
analytically derived [20-22]. A few analytical results have been 
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reported for single [7] and double infections [23]. Farkas & Hinow 
[23] analytically examined the case of hosts infected with only a 
single bacterial strain when two strains were present in the field. 
However, the symbiotic bacterial dynamics of an arbitrary' number 
of strains remains poorly understood [7,18], which prompts the 
questions: "how do the parameters determine the dynamics of 
infections and how many strains can infect a host population"? 
Experimental studies have reported multiple infections of 
Wolbachia in parasitic ants (nine strains) [24], byturid beetles (five 
strains) [25], and ambrosia beetles (five) [26]. However, whether a 
host population can be infected with a larger number of species is 
not known. Thus, the generalization of bacterial numbers is a big 
issue in endosymbiont studies. 

In addition, the modeling of endosymbiotic dynamics in the 
context of mating system has not been extensively studied. CI is a 
phenomenon that occurs via mating, thus it is considered that the 
mating system should affect the dynamics of the endosymbionts. 
Mating systems can be categorized roughly as random mating and 
inbreeding. Random mating or panmictic mating is mating 
between any individuals in the population. In contrast, inbreeding 
or sib mating, is mating between a brother and a sister. Therefore, 
CI should be less effective under inbreeding because the mating 
partner should have the same infection status. Strong inbreeding 
has been observed in some insects, such as fig wasps and ambrosia 
beeties, which are known to have Wolbachia infections [26,27]. 
Some theoretical studies have addressed the dynamics under 
inbreeding [28-30]. However, most of these models were 
examined numerically rather than analytically. In addition, the 
scope of these models was limited in the case of single or double 
infections. 

This study aimed to develop an analytically solvable model of 
the infection dynamics of Cl-inducing bacteria in sib/ panmictic 
mating systems. We derived the analytical equilibrium of single 
and double infections, and generalized arbitrary- numbers of strain 
infections using three parameters (vertical transmission rate, CI 
levels and the level of inbreeding). Our results show that bacterial 
infection thresholds can be simply expressed as a function of the 
ratio between the failure of vertical transmission and the CI level. 
In addition, we derived the msiximum number of bacterial strains 
that are capable of infecting a host population. 

Methods 

Model 

We consider a host population infected with endosymbiotic 
strains. The integer i ( = 0,1, . . . ,2^ — 1) represents the state of 
each host individual, infected or uninfected with bacteria, where i 
is written in base 2 with N bits and each bit represents the state of 
infection with bacteria, i.e., 0 and 1 indicate absence and presence, 
respectively. For example, when i = 5 ( = 101 in base 2), the host is 
infected with the first and third bacterial strains but uninfected 
with the second strain. For convenience, we introduce a binary 
function f{i\,i2), which indicates whether the set of infected 
bacterial strains i\ is included in that of ij. If the set of j'l is a subset 
of that of 12, f(ii,h)=i, otherwise /(/i,!2) = 0. For example, 
/(3,7) = 1 because the set of strains 12 (=111 in base 2) contains 
any strains of i'l ( = 011 in base 2). In addition, we define n{i) as the 
number of strains carried by an individual of type /. 

We develop an difference equation to express the dynamics of 
the frequencies of hosts infected with multiple bacterial strains. Let 
Qi be the frequency of individuals infected with the bacterial set 
For simplicity, we assume that all of the bacteria have identical 
parameters in terms of the vertical transmission rate and CI 
strength among strains. 



Next, we consider the formulation of CI among strains. CI 
results in a decreased hatching rate by a factor 1 — z; therefore, the 
fitness of individual mating incompatibly is reduced to 1 — z. The 

fitness of mating between a male with bacterial strains of i and a 
female with bacterial strains of J is generally described as 

/OV)+(i-z)[i-/ov)]. 

A bacterial strain can be transferred by vertical transmission 
from a mother to offspring at egg stage with a probability of 1 — ^u, 
whereas it fails to be transferred with a probability of fi. In the 
field, measurements of n are typically low (<0.05). An individual 
host infected with the bacterial set j will produce eggs with with a 
probabUity of (1 when f{ij)=\. Note tiiat our 
model assumes a diplodiploid sex determination system in the host 
insects, but the effects of CI can be quite different between 
diplodiploidy and haplodiploidy [31]. 

The occurrence frequency of CI depends on the mating pattern. 
Under random mating, CI occurs depending on the distribution of 
qi. The average fitness of a female with the bacterial set j is 

<lk{f{kj) + {\—z){\—f{kjy)}. In contrast, under inbreed- 

jt=0 

ing, CI depends on ^. Mating partner tends to have the same 
infection status as themselves because the pair are possibly 
produced by the same mother. CI is less frequent under inbreeding 
than random mating because it occurs only under following 
conditions: (i) a mother lacks the bacterial strain(s) because of 
vertical transmission failure and (ii) a father carries the strain(s) 
because of their successful transmission from a common grand- 
mother. It is difficult to formulate the exact effect of CI under 
inbreeding because the frequency in a previous generation is 
needed. For simplicity, we hereafter assume that ji is sufficiently 
small to regard the average fitness of a female with the bacterial set 
7 as 1. Let q' be the frequency of infected individuals in the next 
generation and p the probability of random mating. The 
difference equation of qi is described as follows: 

2^-1 

H= E ?/('V)(i-/*)"*'V"^'-"® 
7=0 

r 1 *" 

(l-rt+f E «WM+(l-zXl-/M» . 

where iv is a normalization factor that maintains the sum of the 
updated host frequencies as unity. Eq.(l) can be rewritten as 
follows: 

2^-1 

7=0 

(2) 

\-pz{-\+ IkfikJ)} . 

k=0 

Eq.(2) is symmetrical for p and z, which means that the effect of 
p on the infection dynamics is equivalent to that of z. In addition, 
Eq.(2) reveals that the dynamics can be described using only two 
parameters, p and pz. If two systems have the same value of pz, 
then the dynamics of these systems shotild be identical. If we use 
z=pz as the effective CI strength, then z should correspond to z in 
the systems with the equivalent dynamics under complete random 
mating (p=l). As a result, Eq.(l) can be rearranged to a more 
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simple form using z as follows: 



H= E ?/OV)(i-/*)""V"^'-"® 

2^-1 

l-z ^ ft[l-/(M] 
4=0 



(3) 



Eq.(3) is identical to the equation developed by Frank [18], but z 
is substituted by z. Frank suggested that if all the bacterial strains 
share the common parameters n and z, the polymorphic 
equilibrium solutions should be symmetric for the frequencies, 
i.e., the frequency of the hosts with the same number of strains 
should be identical. For example, if ^no>0, then 
<l\ 10 = ^loi = ?oi 1 shotdd be satisfied. Frank performed a reduction 
of Eq.(3) as follows: 



k = 0 



(4) 



where Qi is the frequency of a type that carries i different strains, 
in which there are 1^-1 different combinations. We analyze 

Eq.(4) hereafter. 

Frank also proved that a polymorphic equilibrium, where both 
infected and uninfected individuals exist, can be stably maintained 
only if there are host individuals infected with all of the bacterial 
strains in the population [14]. In other words, if the largest n(i) in 
the population is less than N, polymorphism shotdd be lost and 
some bacterial strains wiU eventually fail to achieve infection until 
stable maintenance^ is possible (Figure I). Frank stated that the 
prc'sencc' or absence of individuals infected with all bacterial strains 
in a population determines the fate of the population. Therefore, 
we focus on the dynamics of the frequency of hosts infected with 
all of the bacterial strains in an Af-strain population, gjv- 

Next, we analytically derive the equilibrium for the system 
[Eq.(4)]. We will first focus on simple cases with one or two 
bacterial strains in the population before generalizing the results to 
an arbitrary number of strains. 

Results 

Single infection 

When = 1, the system [Eq.(4)] can be represented as follows: 



(5) 



where w=Qi + Qo{l—zQi) and the total frequency Qo + Qi 
equals unity. 

Because 2i ~ 2l = f ~ 2? + 2l + ^) > there are a maxi- 
w \ zJ 

mum of three fixed points in the system: (2o'2i) = 

(1 0) (^1+^^1-45 i-VT^ 



2 2 

'l-Vl-4a l+Vl-4a 



and 



where c(. = fi/z. The phase portrait 




Figure 1. Schematic showing bacterial infection flows, where 

N = 3. The host population comprises a set of triple infections, three 
sets of double infections, three sets of single infections and a 
noninfection set. The arrows represent vertical transmission failure. All 
sets can exist in the population (enclosed by the solid line) only when 
(/111 >0. By contrast, when qm =0, there are no stable states where the 
frequencies of more than one set of double infection has a positve 
value (enclosed by the dashed line). In that case, the host population 
eventually approaches a stable state that include two bacterial strains 
(encolosed by the dotted-dashed line). The proof was presented in [18]. 
doi:10.1371/journal.pone.0094900.g001 



of the system is shown in Figure 2a. The second and third fixed 

points can exist only when 0 < o( < aj;'' = because a saddle-node 

bifurcation occurs at a = afp. Thus, we can obtain the nontrivial 
fixed points except (1,0) only for 0<a<^. (1,0) and 



stable 



and 



is unstable (Figure 2a). The unstable 



2 ' 2 

l + Vl-4a l-v'l-4a^ 
2 ' 2^ 

fixed point shows the basin boundary. If the initial value of Qi is 
below the unstable fixed point, then Qi decreases asymptotically to 
0, i.e. the population loses the bacterial infection. In contrast, if Qi 
exceeds the boundary due to random drift or other accidental 

i+yr^^ 

events, Q\ can converge towards r . Note that because 

l±4Ei5 

librium, the frequency of infected individuals is certainly more 

than half 

Double infection 

In the case ot N = 2, the system [Eq.(4)] can be represented as 
follows: 



WQ2 



vQi =Q2{l-n)n+Qi{l-n)[l-~z{Qi + Q2 



(6) 



where w= Q2 + 2Qi[l -z{Qi + Q2)] + Qo[\ -z{Q2+2Qi)] and 
the total frequency Qo + 2Qi + Q2 equals tmity. This system has 
five fixed points as follows: (2o'2i>22) = (I'O'O)' 



l-H-^yi-&oi + 2n+n^ l+n+^l-8(X + 2n + n'^ 



.0 , 



l-H+^yi-&a + 2^+H^ l+/i- V^r^8a + 2/I+7? 



,0 , 
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)i - (■ ' r ' ^ \ ^ \ 04 • ^ , ' , ' ^ ^ \ 04 ' T ' — , ' ^ ^ ^ > 

0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 



^5 



Figure 2. (A) Phase portrait where A' = l, /i = 0.1 and z = 0.5. When the initial value of Q\ is less than the unstable equilibrium (open circle), Q\ 
decreases toward zero (extinction of endosymbionts). Otherwise, Q\ moves to a nonzero stable equilibrium point (closed circle on the right side). (B- 
D) Vector fields where N = 1 with /( = 0.05 and z = 0.38, 0.24, and 0.20, respectively. The horizontal and vertical axes represents Qi and Q\, 
respectively. Each arrow represents the difference in Qi and Q\ between a generation, AQ = (Q2 — Q2,Q\ — Qi)- The color of an arrow indicates the 
magnitude of the vector. Solid and open circles indicate stable and unstable equilibrium points. The shaded areas depict the basin of attraction for 
the extinction of endosymbionts. (E-G) Vector fields where ^ = 5 with n = 0.05 and z = 0.75, 0.50 and 0.25, respectively. Each arrow represents the 
difference of Qs and Q4 between a generation, AQ = {Qs —Qs,(^4 — Qi)- 
doi:1 0.1 371 /journal.pone.0094900.g002 



\-2a+ J\-(A+ iCfa-^a^ \-2aL-J\-{A+ii)a-^a?-^, 
^ >0!> I and 



1^ 

2 2 

To examine the stability of the fixed points, we performed linear 

stability analysis for each fixed point. At ( 1 ,0,0), the eigenvalues of 

Jacobian are 1— /^i and (1— /()^. Because the absolute values of 

these eigenvalues are less than 1, the state where there are no 

infected insects are locally asymptotically stable. We numerically 

calculated eigenvalues for the other fixed points for any sets of ji 



and z, and concluded that only 



1 -2o:- ^1 -(4 + /Oo:-4a2 



a. 



l-2o£+ v'l-(4 + A')o<-4a2 



was a stable fixed point and the 



others were unstable. This stable frxed point exists 

■rn (2) , (2) -(4 + /i) + \//f2 + 8^ + 32 

if 0 < a < where a' = — ^ — — • When a. = 

8 



1- 2a +J 1 - (4 + /.i) 0: - Aa^ \-2a- J\-(A + n)a.-Aci?^. 
,0:, I and 



\-2a-J\-{A + ^i)a-A(y? \-2a^J\- (4 + /i)a-4a2 

^ ^ 

2 2 
coalesce due to saddle-node bifurcation as well as when = 1 . 
Because a*,^' < a^P for any value of j.i > 0, the range in which 



Q' 



1.0- 
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Figure 3. Comparison of Ql, in the numerical (colored points) 
and analytical (black curves) solutions where A = 1, 2, 5 and 1 0. 

The parameters used in the simulations were selected from the ranges 
of 0</i<0.05 and 0<z<1.0. The convergence condition was 

<10^"'. The analytically derived equilibria are 

l + v/l~4a. 



shown by a solid line for ^=1 {Q\ 

for N = 2 [Eq.(7)], and dotted lines for A' = 5 and 10 [Eq.(IO)] 
doi:10.1371/journal.pone.0094900.g003 



), a dash-dotted line 
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completely infected host can exist at N = 2 is narrower than that 
when A' = 1 . 

As stated iii the previous section, we assume that is a relatively 
small parameter. By substituting zero for j.i as a zero-order 
approximation, at the stable fixed point and the critical value 
of a^^' are given as follows: 



Q*2- 



l-2a + \/l-4a-4o:2 



■l + \/2 



+ 0(!i). 



(7) 



(8) 



The numerical simulation confirmed that this approximation is 
adequate and that ji had less efifcet on Q2 and a*?' than a. (Figure 3). 



the order of frxed points is 



N-\ N-\ 
=(=-^. When 

Q*^ <Q*^ <Q*^ and the stable fixed point is 2*"- However, 
we did not observe a transcritical bifurcation for any parameter set 
of A^, z and /i. Instead, a discrete jump into 0 was observed before 
Q*f^ reached to 0, which implied that a saddle-node bifurcation 
also occured at A^>3 (Figure 3). The reduced system did not 
rephcate the disappearance of the stable frxed point; however, the 
analytical solution of Q*f^ agreed with the numerical results until 
the discrete jump occured. In addition. The numerical simulation 
confirmed that when the bacterial strain number was 
sufficiently high, a discrete jump occured around Qm = 0. 
Therefore, with a large A' limit, the value of a at which the 
transcritical bifurcation occurs is equal to a at the saddle-node 
bifurcation. Thus, we obtained the critical aj.^* for a sufficiendy 
large A^ as follows: 



A'^-strain infection 

Finally, we calculated the frxed points of Eq.(4) for arbitrary A^. 
We assumed that Q„ x 0, for 0 < « < A^ — 2 in an equilibrium state. 
Numerical simulations supported this assumption (Figure SI). For 
example, when Af=10, 5 = 0.4 and /( = 0.025, 95.5% of the 
population possessed A' or A^ — 1 bacterial strains. The dominance 
of highly infected hosts was enhanced by the strain number A'^. 
Thus, this assumption is more plausible when the system has a 
greater A' and a lower a. The reduced system is described as 
follows: 



JN). 



(9) 



where Qn + NQn-i^I, wx Qm + NQf^^i[l -~z(Qm + (N -I) 
Qn-i)]- Because this is a one-dimensional system, the fixed points 



2-Af- 



N[- 4 + Nz+ 4(1 -{.if 



can be found easily, 2]^ = 0, 

If l-i was regarded as a relatively small parameter, 
{1-f.if x\-Nii. Then, 



2^=0, 



2-N±N\/l-4a 



(10) 



which shows that the frequency of completely infected hosts 
depends only on a when /( is sufficiently small. The difference in 
giv between one generation, AQ^j, is represented as follows: 



AQn = Q'n -Qa 



wN 



(Qn - Q°)iQN - e+ )iQN -Q-),W 



where Q*° = 0, Q*^ = ^ N±NV\ 4a ^^^^ 

a < , the order of the three fixed points is 

Q*^ <Q*^ <Q*^ and Q*^ is a stable equilibrium. 

Figure 3 compares Q*^ at various a with the results calculated 
by numerical simulations for A' = 5 and 1 0. The approximated 
solution of Q*"*" was generally consistent with the numerical 
results. In particular, the agreement was better when a was 
sufficiently small. 

Eq.(ll) indicates that an exchange of stability, so-called 
transcritical bifurcation, should occur between 2*+ and Q*^ at 



A^-l 



(12) 



We compared the analytically derived a.c with numerical 
simulations. Figure 4 presents the numerically calculated Q*f^ at 
various /.i and z when A^= 1, 2 and 5. The lines in each figure 

indicate the analytically derived critical value of (X, — = a;. ' = -, 




z 



Figure 4. Distribution of the values of obtained by 
numerical simulations where A' =1,2 and 5. The darl< highlights 
indicate the values of gw- ^'^^ Wries in each figure represents the 
analytically derived critical value of a. 
doi:10.1371/journal.pone.0094900.g004 
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a 0.25- ■ 




-1 — I — I — I — I — I — I — I — I — I — I — r 
2 4 6 8 10 12 

Figure 5. Threshold of a. that allows bacterial infection. The 

numerical simulation confirmed that a host population where a is in the 
shaded region could maintain infections of N bacterial strains. The 
black crosses are the analytically derived thresholds. The horizontal 
dashed line is drawn at a = 0.25. 
doi:1 0.1 371/journal.pone.0094900.g005 



v(2)- 



■1 + 72 



and d 



(5). 



-. The lines successfully delineate the 



2 25 
parameters where bacterial infections can be maintained. The 

consistency is reduced when \i is greater because a[^' and 
include 0(/() errors. In particular, the ranges of z and n that allow 
infections with a complete set of bacteria are distinct. The value of 
z ranges from 0 to 1. In contrast, when fi>!Xc, Qn cannot have a 
positive value, irrespective of the value of z. Figure 5 confirms the 
consistency between the numerical results and Eq.(12). The shaded 
area derived by simulation indicates the values of a where A'^ 
bacteria can infect a population. On the basis of the analytical 
results, we plotted the critical value of a where A'^=l and 2, 

"c'^^^) — ^'2^ ' ^ — 3) we plotted the critical a 
derived from Eq.(12). For larger values of A^, the analytical results 
approached the threshold derived by numerical simulation. 

The inverse function of Eq.(12) defines the maximum number of 
bacterial strains that can infect a population, Nc, which is given as 
follows: 



Nc 



I + v'l-4a_ 1 
2a ~ a 



-1. 



(13) 



Nc is approximately proportional to a ' . Therefore, it has a larger 
value with a lower value of a. 

Discussion 

In this study, we demonstrated the analytic equilibrium 
solutions of the frequencies of individuals infected with all the 
bacterial strains in a population with single, double and arbitrary 

infections. Our results provide of a qualitative insight into the 
symbiotant bacterial dynamics, in constrast to recentiy developed 
models that quantitatively simulate specific experimental results 
[4,32]. 

Our model used a parameter p, the level of inbreeding, from the 
perspective of host behavior. We revealed that p has a completely 
negative efiect on CI. Endosymbionts manipulates host reproduc- 
tion to maintain their own infections in host populations. In 
contrast, host insects can reduce the prevalence of endosymbiont 
infections by inbreeding. This result is consistent with previous 



numerical simulations of inbreeding [28]. Engelstadter et al. 
claimed that the infection rate decreases with increasingly inbred 
hosts because uninfected females increasingly mate with uninfect- 
ed males, which leads to fewer incompatible matings. Our analytic 
calculation determined a linear relationship z=pz, which clearly 
indicates the negative effect of inbreeding on the CI strength. 
However, this simple relationship does not appear to hold in some 
cases of stochastic models. Genetic drift in stochastic island model 
might have effects on inbreeding [29]. Our result is also consistent 
with the previous analytical result reported by Dannowski et al. 
[30]. This model focused on double infections by mak'-killing 
bacterial strains and analytically showed that higher level of host 
inbreeding leads to lower frequency of bacterial infections. 
Although their result was limited in the case of double infections, 
our result is applicable to arbitrary number of Cl-inducing 
bacterial strains. 

The analytically derived erjuilibrium depends mainly on one 
parameter, a. a is a novel index that we defined in this study, 
which is the ratio of the probability of vertical transmission failure 
(fi) relative to the effective CI level (z). Previous studies have shown 
that these parameters are involved with the infection dynamics. 
Furthermore we demonstrated the simple relationship between 
these parameters, i.e., the bacterial infection dynamics are 
determined simply by a balance of both. The equilibrium 
frequencies of the fuUy infected host are at a higher rate when a 
is lower (Figure 3). In this case, Cl-inducing bacteria can select two 
alternative strategies to maintain and spread themselves: increas- 
ing the vertical transmission efficiency [16] and/or increasing the 
CI strength. There have been no reports of the differences among 
these strategies. Figure 4 presents the asymmetric relationship 
between fi and z. When fi was higher than approximately 0.2, 
bacterial infections could not be sustained regardless of the value 
of z. Thus, the vertical transmission rate has a threshold value for 
infection. In contrast, bacterial infections occur with any value of z 
provided an appropriate value of p is selected on the basis on z. 
Thus, the vertical transmission rate determines the infection state 
more strongly than CI. 

Eqs(12) and (13) represent the relationship between the 
maximum number of bacterial strains and the bacterial param- 
eters. Considering a host population containing A'o bacterial 
strains, if A'o inc:reases to A^i through the horizontal transfer of a 
bacterial strain, the critical valule of a decreases from aj/^"' to a'^''. 
Provided a<a*^'' is satisfied, the A^i -infecting population will be 
maintained. Otherwise, the infection cannot be maintained. For 
example, a population where a = 0.2 can maintain three strains at 
most because ~ 0.1 8) <o;<o;*''(~ 0.22). If the population 
actually contains the maximum number of strains, the host 
population is "saturated" with bacteria, i.e. other strains will fail to 
invade the host population. The parameters are expected to have 
the same value in our model. Therefore, the frequency of the 
rarest strain necessarily declines to zero. A higher a value than that 
of other strains may be needed for a new strain to invade the 
saturated hosts. 

Our theoretical results suggest that the equilibrium frequency of 
completely infected hosts decreases as N increases (Figure 3). In 
addition, the range of a where bacteria can be maintained 

becomes narrower as N increases (Figure 5). The severe conditions 
for multiple infections are caused by the increased possibility of 
vertical transmission failure the offspring by any strains. X-strains 
that infect a individual can successfully transmit aJl the strains to 
the offspring with a lower probability than X—l strains. 
Consequently, these results support a noble idea of Wolhachia- 
induced speciation. Suppose that a is enough large to maintain 
single infection but not double infection (0.207-0.25), and 
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Figure 6. Comparison between the analytic equilibria for A' = 1 

and 2, and the reported data. The solid and dashed lines represent 
the analytic estimates for ^ = 1 and 2, respectively. The closed circles 
and an open square indicate each single infection and double infection 
dataset, respectively. The values of /i and z, and the references are 
shown in table SI. 

doi:1 0.1 371/journal.pone.0094900.g006 

endosymbiotic strains, A and B strains, infect some insects in two 
diflFerent geographic regions respectively. The infections in both 
regions will spread until a hybrid zone is produced. Then, 
horizontal transfer makes double-infected individuals of A and B 
strains. However, the double infection must be lost because a is too 
low to maintain the both infections. Therefore, the hybrid zone 
wiU be unstable and post-zygotic speciation wiU be occurred 
between A-infected and B-infected individuals. 

To confirm that our results were consistent with experimental 
data, we compared the predicted equilibria with the values of CT 
inducing Wolbachia and Cardinium obtained using natural popula- 
tions or an artificial line (Figure 6, Table SI). No data were 
available for triple or higher multiple infection. Therefore, only 
five single infections and one double infection were used for 
reference. No data were available on the inbreeding frequency of 
the insects studied; hence, p values of all were assumed to be zero, 
i.e. completely random mating. The experimental values were 
close to the predicted fines. Thus, our results can be used to 
estimate the frequency of infected hosts at equilibrium using the 
parameters [p. and z) as well as for estimating the parameters that 
are often difficult to be measured. 

In this study, we aimed to provide a qualitative outiine of the 
dynamics of multiple infections by introducing an analytically 
solvable model. By ignoring quantitative accuracy, our model 
sheds light on the mathematical structure related to the multiple 
infection dynamics. However, we introduced several assumptions 
to develop a solvable model, which should be eliminated to allow 
more quantitative comparisons using experimental data. First, we 
assumed that a common transmission rate or CI level were shared 
among bacterial strains. If the bacterial parameters are heteroge- 
neous, a bacterial strain with a lower a is more likely to be 
maintained than other strains. Indeed, some experiments indicate 
that the CI range is too broad to be regarded as the same value 
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